Phase transitions in a disordered system in and out of equilibrium 
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The equilibrium and non-equilibrium disorder induced phase transitions are compared in the 
random-field Ising model (RFIM). We identify in the demagnetized state (DS) the correct non- 
equilibrium hysteretic counterpart of the T = ground state (GS), and present evidence of univer- 
sality. Numerical simulations in d = 3 indicate that exponents and scaling functions coincide, while 
the location of the critical point differs, as corroborated by exact results for the Bethe lattice. These 
results are of relevance for optimization, and for the generic question of universality in the presence 
of disorder. 

PACS numbers: 



Similarities and differences between equilibrium and 
non-equilibrium states in disordered systems have been 
widely studied both for their conceptual importance and 
because the presence of randomness often provides pro- 
totypical examples of complex optimization problems . 
There are also many applications in the physics of ma- 
terials, where this dichotomy is met, together often with 
concepts such as aging and glassiness. The central issue is 
to understand whether the equilibrium properties of dis- 
ordered systems provide a faithful representation of the 
non-equilibrium states in which the system is likely to be 
found in practice. In optimization terms, the question is 
what is the relation between an approximate solution and 
the optimal one. 

A disordered system can be non-trivial even at zero 
temperature due to the presence of a complex energy 
landscape. The properties of the ground-state (GS) are 
often difficult to determine analytically, and numerical 
evaluation becomes computationally prohibitive for large 
systems in particular and for some problems like spin 
glasses, in general. Non-equilibrium dynamics brings the 
system to the nearest metastable state and then noise or 
an applied field is needed to allow it to explore the en- 
ergy landscape. Typical optimization methods are con- 
structed by providing a suitable perturbation scheme on 
the states of the system. Recently hysteretic optimiza- 
tion was proposed 0] as an alternative to methods that 
use noise, such as simulated annealing 0. Its basis is an 
analogy to a ferromagnetic demagnetization procedure: 
an external oscillating field with decreasing amplitude 
and low frequency is applied to the system, yielding at 
zero field the demagnetized state (DS), which is used as 
a reference state for material characterization 

The ferromagnetic random field Ising model (RFIM) 
has been extensively studied in literature as a paradig- 
matic example of disordered system whose equilib- 
rium and non-equilibrium properties are still tractable, 
though far from trivial. The RFIM is one of the simplest 
systems, where the crucial interplay between quenched 
disorder and exchange interaction gives rise in high 
enough dimensions to a disorder induced phase tran- 
sition. This also affects the dynamics, where a non- 



equilibrium phase transition exists. 

The equilibrium properties of the RFIM are governed 
by the T — scaling Q even at high temperatures. GS 
calculations have elucidated the properties of the phase 
diagram: In d = 1 the RFIM is trivially paramagnetic. 
In d = 2 there is no phase transition but an effective 
ferromagnetic regime for small systems. In d > 3 the GS 
displays an equilibrium phase transition induced by the 
disorder from a low disorder regime where the system 
is ferromagnetic (FM) to a strong disorder one, where 
the system is paramagnetic (PM) 0, 0] • The equilibrium 
critical exponents for random field magnets have been 
measured experimentally in Feo.93Zno.07F2 00] 

Likewise, the non-equilibrium properties of the RFIM 
have been studied by extensive numerical simulations and 
renormalization group calculations 0,0, El- A disorder- 
induced transition is observed in the hysteretic behavior 
for d > 3: At low disorder the loop has a macroscopic 
jump in the magnetization, which disappears at a criti- 
cal value of the disorder, above which the loop is smooth 
on a macroscopic scale. Numerical simulations ^3 an( i 
renormalization group have been used to estimate the 
critical exponents in various dimensions. A disorder in- 
duced non-equilibrium phase transition in the hysteresis 
loop has been studied experimentally in C0-C0O films 
[Hi and Cu-Al-Mn alloys [H . 

The relation between this non-equilibrium transition 
and the PM to FM one in the (equilibrium) GS has been 
debated in the past. Based on the similarity in the nu- 
merical values of the exp onents and on mean-field equa- 
tions, Maritan et al. |l3( argued that the two transitions 
should be universal. Sethna et al. 0] refuted this due 
to the different natures of the two cases: the transition 
in the GS occurs for a zero external field, while the tran- 
sition in the hysteresis loop occurs at the coercive field. 
More recently, the question of the universality of the ex- 
ponents, with respect to the shape of the disorder dis- 
tribution, was discussed in d = 3 simulations, mean-field 
theory, and on the Bethe lattice 0, 0, 0] . 

In this letter, we compare the equilibrium and non- 
equilibrium phase transitions in the RFIM, with evidence 
for universality. The key issue is the identification of a 
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reference non-equilibrium state, instead of focusing on 
the jump in the saturation loop. In the low disorder 
phase, a discontinuous hysteresis loop corresponds to a 
region of the field-magnetization plane not accessible by 
any field history |18l [lSj . In this regime it is not possible 
to demagnetize completely by applying a slowly varying 
AC field. Thus one studies the DS as the state of lower 
(remanent) magnetization resulting from the demagneti- 
zation procedure. This state is uniquely defined, in the 
quasistatic limit, for any given realization of the random 
fields. It has two non-equilibrium phases: FM when the 
main loop has a jump, and PM otherwise. The rema- 
nent magnetization becomes the order parameter of the 
transition. Notice that the DS is defined at H — and 
is therefore the natural non-equilibrium counterpart of 
the GS. This responds to the objection raised in Ref. 0] 
against the possible existence of universality. 

We evaluate in d = 3 the finite-size scaling functions 
both for the equilibrium and non-equilibrium phase tran- 
sition. By rescaling the disorder around the appropriate 
(distinct) critical values, the scaling functions can be col- 
lapsed onto the same curve using the same exponents 
values. However, the location of the critical point differs: 
the transition in the DS occurs at a lower disorder value. 
Thus there is an intermediate region where the GS is fer- 
romagnetic but the DS is paramagnetic. To further cor- 
roborate our findings, we analyze the RFIM on the Bethe 
lattice: we compute the GS magnetization and compare 
it with the remanent magnetization of the DS 0. While 
the exponents are the same in the two cases (coinciding 
with mean-field results) , the Bethe lattice reproduces the 
ordering of the critical points in d = 3. Additional evi- 
dence for universality is obtained by comparing the order 
parameter distribution function at the critical point for 
finite systems. 

In the RFIM, a spin Si = ±1 is assigned to each 
site i = l,...N (here of a cubic lattice in d = 3 or 
a Bethe one with coordination z = 4). The spins 
are coupled to their nearest-neighbors spins by a fer- 
romagnetic interaction of strength J and to the external 
field H. In addition, to each site is associated a ran- 
dom field hi taken from a Gaussian probability density 
p(h) = cxp(—h 2 /2R 2 )/V2TrR, with width R, which mea- 
sures the strength of the disorder. The Hamiltonian thus 
reads 

H = -Y J Js i s j -^{H + h i )s i , (1) 

where Yl/j j\ is restricted to nearest- neighbors pairs. 

The RFIM GS is numerically solvable in a polynomial 
CPU-time with exact combinatorial algorithms. We find 
the GS via the min-cut/max-flow problem of combinato- 
rial optimization, and use the so-called push-relabel vari- 
ant of the preflow algorithm [2(|. Such methods, properly 
implemented, are in general slightly non-linear in their 
performance as a function of the number of spins [lj . 

For the out of equilibrium case, we consider a simple 
relaxation dynamics obtained in the limit T — > of the 



Glauber dynamics 0, 0, : At each time step the spins 
align with the local effective field 

Si = sign( Sj + hi + H) (2) 
i 

until a metastable state is reached. To construct the hys- 
teresis loop, the system is started from a state with all 
the spins down = —1 and then H is ramped slowly 
from H -> -co to H oo. The limit of dH/dt -> 
is taken after the limit T — > 0. In practice, this can be 
conveniently obtained by precise increases of the field, 
to always flip the first unstable spin. To reach the DS, 
the external field is changed through a nested succession 

H = H -» H Y -> H 2 -> H n ... -» 0, with H 2n > 

—H 2n +i > H 2n +2 > 0, and dH = H 2n — H 2n+2 — > 0. 
This provides a perfect demagnetization with a uniquely 
defined DS (dH/dt — ► 0). It being quite expensive com- 
putationally, we instead perform an approximate demag- 
netization using an algorithm discussed in Ref. |18( with 
dH = 1CP 3 . We verified that the states have negligible 
differences with perfectly demagnetized ones. 

The RFIM critical exponents characterizing the disor- 
der induced transition can be defined as usual: The mag- 
netization M = (\m\), with m = J2i S i/N, scales close to 
the transition point as M — Ar@ , where r = (R — R c )/R c 
is the reduced order parameter and A is a non-universal 
constant. The correlation length exponent £ = (Br)~ v , 
where B is another non- universal constant, rules the 
finite-size scaling of the model 

M = AL-V v f (bi}I v (R - R c )/R c ) . (3) 

GS simulations in d = 3 for Gaussian disorder, yield 
Ujy(GS) _ 073j p(GS) _ o.02 and R ( C GS) ~ 2.28 

[HmiHlil. 

The demagnetization process has been exactly solved 
in d = 1 ^3 and on the Bethe lattice (0|> see a ^ so 
|25|). Numerical simulations in d = 3 indicate that the 
DS displays the same critical point as the saturation loop 
|2(j| . The transition point has been obtained numerically 

in d — 3 as Ri DS ^ ~ 2.16 ^(| an d the critical exponents 
have been measured. E.g. Ref. [2(| reports (3(ds) = 0.04± 
0.02 and l/v {D s) = 0.71 ± 0.1. 

The numerical simulations for the GS and DS are done 
for the same disorder realizations for the both cases, for 
cubic lattices of linear sizes L — 10, 20, 40, 80. The results 
are averaged over several realizations of the quenched 
random fields. In both cases, we compute the average 
magnetization as a function of the disorder width. In 
Fig. n we collapse the two sets of data into a single curve, 
using two different values for R c (i.e. Ri GS ^ — 2.28 and 
r( ds ^ = 2.16) but the same values for the exponents (i.e. 
1/v = 0.73 and j3 = 0.03). The best value for the ratio of 
the non-universal constant is found to be Ads/Aqs ~ 1 
and B DS /B GS = 0.68 ± 0.02. 

To provide another viewpoint and corroborate our 
claims, we compare the GS and DS on the Bethe lat- 
tice where analytical expressions can be found exactly. 
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FIG. 1: Numerical results in d = : The magnetization can be 
collapsed using R c = 2.28 (GS) and R c = 2.16 (DS), l/v = 
0.73 and /3 = 0.03. The scaling curve is the same for DS 
and GS indicating universal behavior. The values for the 
ratios of the non- universal constants are Ads/Ags = 1 and 
B D s/Bgs = 0.68. 



FIG. 2: The magnetization of the GS and the DS computed 
exactly on the Bethe lattice with z — 4 in the thermodynamic 
limit, showing the ordering of the critical point (see inset). 
When the data are plotted against the reduced parameter 
(R c — R)/R c the curves superimpose. The result implies that 
for the Bethe lattice Acs = Ads- 



The RFIM displays also on the Bethe lattice, for a large 
enough coordination number z, both an equilibrium and 
a non-equilibrium disorder induced phase transition [l^ . 
To compare the GS and the DS around the respective 
transitions, we take directly the thermodynamic limit, 
using for the DS the results of Ref. We have ob- 

tained the GS magnetization following Refs. 27, 28] as 



M 



/OO p OO z 

dhp(h) / Y[dx k Woo(xk)(s ). (4) 
- OO J — OO 



Here W cx> (x) is fixed-point probability distribution for 
the quantity x n = -j\n(Z+ /Z~), where are defined 
as the partition functions of a branch of generation n 
with a fixed up (down) spin so at the central site [13, EU, 
and is given by an implicit integral equation. The fixed 
point equation is solved by numerical integration, and 
the magnetization is computed for different values of R 
using Eq. fljl. for T = and z = 4. In Fig. we show a 
comparison between the magnetization in the GS and in 
the DS ([13). 

As in d — 3 simulations, the transitions occur at two 
different locations (see the inset of Fig. [2J, for z = 4 
R { C DS} = 1.781258... [3 and R { C GS) ~ 1.8375, with the 
mean- field exponent ((3 = 1/2). When plotted against 
(R — R c )/R c the two curves superimpose close to the 
critical point. This indicates that, though not required 
by universality, in the Bethe lattice Acs = ^ns, as also 
found in d = 3. To investigate possible finite size scaling 
we have performed numerical simulations in the Bethe 
lattice, following the method of Ref. 29] . Collapsing the 
order parameter curve as in d — 3, using a scaling form 
similar to Eq. @, does not appear to be possible in the 



Bethe lattice, because the scaling region is very narrow. 
Thus to test finite size scaling, we have computed the dis- 
tribution of the magnetization m at the respective critical 
point, i?c £l ' S ' ) and r{ GS,> for different lattice sizes N. The 
distributions can all be collapsed into the same curve (see 
Fig. E|, using the form P(\m\) = f(\m\/M)/M. 
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FIG. 3: The distributions of the magnetization in the DS and 
the GS at their respective critical points on the Bethe lattice, 
obtained numerically for different lattice sizes N, can be all 
collapsed together. 

To conclude we provide evidence about the universality 
of the RFIM with Gaussian distributed disorder in and 
out of equilibrium. The key point is the identification 
of the correct order parameter for the non-equilibrium 
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transition. This quantity, the remanent magnetization of 
the DS, is the natural counterpart to the magnetization 
of the GS in the equilibrium case, in particular at a zero 
external field. Our results are based on a detailed nu- 
merical analysis in d — 3 and on an exact solution on 
the Bethe lattice. It would be interesting to confirm this 
conclusion by more complex measures, beside the order 
parameter, such as the scaling of the domain wall stiff- 
ness [24j . Regardless of the question of universality, the 
most intriguing point from our analysis is that the two 
transitions appear at different critical values of the disor- 
der strength: the DS ferromagnetic phase is the first to 
disappear as the disorder is increased (Ri DS ^ < Ri° S ^). 
The interpretation of the ordering is simply that the GS 
is as correlated as a state can be in the RFIM, due to 
the global optimization. Thus eg. as R is decreased it is 
natural that the FM correlations appear in the GS first. 

Our results have important consequences on the use 
of the demagnetization as an optimization tool: the dif- 
ference in the location of the critical points implies that 

for Ri DS ^ < R < r{ GS \ the DS (paramagnetic) is dras- 
tically different from the GS (ferromagnetic), suggesting 
that in that region hysteretic optimization is likely to 
fail. Moreover, to compare the GS and DS in a system, 
one expects to achieve the closest resemblance in this 
regime if the correlation length is the same; ie. due to 



the difference of the critical points at two separate val- 
ues Ri < R^ S \ respectively - or as well two values 
of the effective coupling constant J in the Hamiltonian. 

In addition to the ferromagnetic RFIM model, one 
can speculate about other systems where two disorder 
induced phase transitions exist. Numerical simulations 
and analytical results have shown that a disorder induced 
transition in the hysteresis loop can be observed in the 
random bond Ising model |3(| . in the random field O(N) 
model [3 11], in the random anisotropy model [32( and in 
the random Blume-Emery-Griffith model |30|. All these 
systems also show a transition in equilibrium and it would 
be interesting to compare their DS and GS. Another ex- 
ample would be the study of an interface in quenched dis- 
order, where many results are known for the roughness 
exponent in and out of equilibrium (i.e. at the depinning 
threshold), and the results typically differ [jj. It would 
be interesting to measure the roughness of an interface 
after a demagnetization cycle (i.e. after the field driving 
the interface is cycled with decreasing amplitude), and 
compare its properties with those of the ground state 
interface. Finally, there is the issue of energetics of ex- 
citations in the respective ensembles: the universality of 
exponents and scaling functions would seem to imply that 
these also scale similarly. 



[16] 



M. Alava, P. Duxbury, C. Moukarzel, and H. Rieger, 
in Phase transitions and critical phenomena, Vol 18, [17 
edited by C. Domb and J. Lebowitz (Academic Press, 
San Diego, 2001). [18 
G. Zarand, F. Pazmandi, K. F. Pal, and G. T. Zimanyi, 
Phys. Rev. Lett. 89, 150201 (2002). [19 
S. Kirkpatrick, CD. Gelatt, and M.P. Vecchi, Science 
220, 671 (1983). [20 
G. Bertotti, Hysteresis in Magnetism (Academic Press, [21 
San Diego, 1998). [22 
T. Nattermann, in Spin Glasses and Random Fields 
edited by A.P. Young (World Scientific, Singapore, 1997). [23 
Z. Slanic, D. P. Belanger, and J. A. Fernandez-Baca 
Phys. Rev. Lett. 82, 426 (1999). [24 
F. Ye et aL, Phys. Rev. Lett. 89, 157202 (2002). 
J. P. Sethna et al, Phys. Rev. Lett. 70, 3347 (1993). [25 
K. Dahmen and J. P. Sethna, Phys. Rev. B 53, 14872 [26 
(1996). 

O. Perkovic, K. A. Dahmen, and J. P. Sethna, Phys. Rev. [27 
B 59, 6106 (1999). [28 
A. Berger, A. Inomata, J. S. Jiang, J. E. Pearson, and 
S. D. Bader, Phys. Rev. Lett. 85, 4176 (2000). [29 
J. Marcos, et. al. Phys. Rev. B 67, 224406 (2003). 
A. Maritan, M. Cieplak, M. R. Swift, and J. R. Banavar, [30 
Phys. Rev. Lett. 72, 946 (1994) [31 
J. P. Sethna et al, Phys. Rev. Lett. 72, 947 (1994) 
J.-C. Angls d'Auriac and N. Sourlas Europhys. Lett. 39, [32] 
473 (1997). 

P. M. Duxbury and J. M. Meinke Phys. Rev. E 64, 



036112 (2001). 

R. Dobrin, J. M. Meinke, and P. M. Duxbury, J. of Phys. 
A: Math. Gen. 35, L247 (2002). 

L. Dante, G. Durin, A. Magni, and S. Zapperi, Phys. 
Rev. B 65, 144441 (2002). 

F. Colaiori, A. Gabrielli, and S. Zapperi, Phys. Rev. B 
65, 224404 (2002). 

E. Seppala, PhD thesis, HUT, Helsinki, 2001, 

A. T. Ogielski, Phys. Rev. Lett. 57, 1251 (1986). 

A. K. Hartmann and U. Nowak, Eur. Phys. J. B 7, 105 

(1999). 

A. K. Hartmann and A. P. Young Phys. Rev. B 64, 
214419 (2001). 

A. A. Middleton and D. S. Fisher, Phys. Rev. B 65, 
134411 (2002). 

P. Shukla, Phys. Rev. E 63, 027102 (2001). 

J. H. Carpenter and K. A. Dahmen, Phys. Rev. B 67, 

020412 (2003). 

R. Bruinsma, Phys. Rev. B 30, 289 (1984). 

M. R. Swift, A. Maritan, M. Cieplak, and J. R. Banavar, 

J. Phys. A 27, 1525 (1994). 

D. Dhar, P. Shukla, and J. P. Sethna, J. Phys. A 30, 
5259 (1997). 

E. Vives and A. Planes, Phys. Rev. B 50, 3839 (1994). 
R. da Silveira and M. Kardar, Phys. Rev. E 59, 1355 
(1999). 

E. Vives and A. Planes, Phys. Rev. B 63, 134431 (2001). 



